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Abstract 

We present a computational study using Molecular Dynamics of the devel- 
opment of an accumulation of grains in one side of a two dimensional U-tube 
under vertical vibrations. Studying the evolution of the height difference 
between the centers of mass of the branches of the tube, we found that it 
reaches a saturation value after vibrating for some time. We obtain that this 
saturation value is the same if the simulation start with the arms leveled or 
with a large initial height difference. We explore the effect of the width of 
the tube, the density of the grains and the coefficient of restitution between 
the grains and the wall on the value of the saturation. We obtain a value 
of the width of the tube for which the saturation value reaches a maximum, 
and show that the transport of grains between arms is favored for low grain 
densities and high grain- wall restitution coefficient. 
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1. Introduction 

Vertical vibrations can induce a collective motion of grains inside a U- 
shaped tube or a partitioned rectangular container. The main consequence of 
this collective motion is the accumulation of the granular material in one side 
of the tube, that can depend on several parameters like the shape of the con- 
tainer, the properties of the granular material, the frequency and amplitude 



Preprint submitted to Advanced Powder Technology 



March 8, 2013 



of oscillation and the presence of any interstitial fluid. When the effect of 
an interstitial fluid between grains is relevant, the height difference between 
the limbs of a vibrated U-tube tends to grow exponentially. This behaviour 
has been studied experimentally for small grains (~ 200 fim diameter) in air 
[3, 0, and theoretically, computationally and experimentally for coarser 
grains (~ 800 fim diameter) in water [ij. 



Akiyama et al. [5] instead of using an U-tube, studied experimentally a 
rectangular container partitioned by a vertical wall on its center, with a small 
gap below for communication between both sides of the container. They ob- 
served that vibration on an initially leveled partitioned container, induced 
the accumulation of grains on one side of the container, in a way that the 
level difference between the partitions grew up to a certain height. They 
showed that the maximum height difference achieved was larger when the 
experiment was performed at atmospheric pressure, but shrank to a value 
of ~ 20 grain diameters, when the air pressure was reduced below 50 Torr. 
They also presented a simple snapshot of a two-dimensional molecular dy- 
namics simulation of soft spheres, showing a maximum height difference in 
accordance with the experimental result at reduced air pressure, however no 
further information is given regarding the simulations. 

Ohtsuki et al. [6;], studied a partitioned container like Akiyama et al. 
but with unequal partition widths. They also showed simulations to com- 
plement their experimental results, using two-dimensional polydisperse hard 
discs. In their simulations, the accumulation of grains always occurs towards 
the widest partition, while it can occur towards any side in symmetrically 
partitioned containers or U-tubes. They focused their computational study 
on how the kinetic friction between the walls and grains affects the transport 
of grains from one partition to another. 

Until now, previous research has identified two key factors affecting the 
development of a height difference in partitioned containers or U-tubes: En- 
ergy loss due the grain-wall interaction and the presence of an interstitial 
fluid between grains. With the help of Molecular Dynamics simulation 0, S] 
we decided to study the problem in the absence of an interstitial fluid, and 
in a simple two dimensional case, to characterize the problem in a controlled 
basic configuration. We studied the evolution of the height difference of the 
center of mass of the U-tube arms Ay, and also the dependence of the sat- 
uration value Ay sat , as a function of the tube width, the mass of the grains 
and the coefficient of restitution between the walls and the grains. 
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Figure 1: Model of the granular material inside a U-tube under vibration. 

2. Computational model 

Our model consist of N monosized discs with mass m and diameter d, 
partially filling a thin U shaped tube with vertical branches of width d t , 
height L y and external width L x (see Fig. [1]). The tube is vertically vibrated 
applying an acceleration of the form a y w 2 sen(ut) , where a y is the amplitude 
of vibration, w is the angular frequency and t is the time. 

In our system, two grains i and j with positions {r*j, fj} and translational 
velocities {vi, Vj}, will suffer an inelastic collision only when the separation 
between their centers of mass is lower than their effective diameter = 
\fi — fj\ < (di + dj)/2. During the collision, the ith grain will feel a restorative 
and dissipative contact force due to the jth grain, with normal component 
(through the line connecting their centers of mass, n direction) and tangential 
(through a line perpendicular to the line connecting their centers of mass, 
passing through the contact point, s direction) of the form: 



Fij = {k n (dij - rij) - 7„m r u"}n - — ^ — - {min{~j s m r \vlA, fi\Fij ■ n\)}s, (1) 

\ v ij ' s \ 
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where ?% = v { - v j} h = r^/r^, v% = ■ n, s = and vfj = 

Vij — v^n. The first term on the right represents the elastic grain-grain 
interaction acting during contact, the second term describes viscous forces 
for the normal deformation proportional to the collision velocities, and the 
last term accounts for the tangential friction between grains, where we take 
the minimum between viscous and Coulomb's friction. Thompson and Grest 
j9[ found that rotations were irrelevant in two-dimensional modeling of a 
monosized granular material in the presence of gravity, therefore as first 
approach we do not include rotations in our model. 



Tabic 1: Parameters used in simulations in units of m, g and d. 

Number of grains 610 — 2210 

Grain diameter 1 

Grain mass 0.25 — 2.5 

Grain density 0.32 p — 3.16 po 

Elastic constant k n 7989 k - 199992 k 

Normal dissipation coefficient 7„ 7.98 7 — 1974 70 

Tangential dissipation coefficient j s 7.98 70 — 1974 70 

Friction coefficient p 0.5 

Grain-wall restitution coefficient e w 0.05 — 1 

Tube width d t 5.5 d — 16.67 c? 

Tube height L y 120.5 d 

Base width L x 23.1 d - 45.44 d 

Vibration amplitude a y 6.83 d 

Vibration frequency w 0.49 wo 

Time step At 10~ 4 1 



The parameters k n , 7„ iS and fi, are the elastic constant, the coefficients 
of viscous friction and the kinetic coefficient of solid friction and m r = 
miTrij/(mi + rrij) is the reduced mass. The corresponding contact force on 
the jth grain is given by the third Newton's Law = —F^. Under the 
gravitational field g, the translational accelerations of the grains are deter- 
mined by the second Newton's law, in terms of the total force acting on 
grain i: F\ ot = rrtig + F^. To save computational time during the calcu- 
lation of the interaction force between grains, we used the cell algorithm jsj. 
Part of the energy lost during a collision is quantified through the restitution 
coefficient e n , for the normal deformation, given by e n = e~ lntco ^ 2 , where 
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Figure 2: Temporal evolution of TV = 1700 grains in a U-tube, filled up to L y /2 with 
L x = 35.7 d and d t = lid. Starting with the arms leveled at time zero, the height 
difference Ay between the centers of mass of each arm grows to a saturating value. The 
grain accumulation can happen in any arms of the tube, as it happens in the experiments 
at reduced air pressure reported by Akiyama et al [5] . 



tcoi = K/(k n /m — 7^/4) 1 / 2 , is the collision time, and the elastic constant k n , 
is chosen such that two grains do not overlap more than 1% under the ac- 
tion of gravity. Taking this in consideration, in order to guarantee numerical 
stability our time step was taken as At ~ t co ;/100. 

The U-tube has rough walls, such that the normal and parallel compo- 
nents of the pre-collision Vi and post-collision grain velocities are given 
hy v' i± = — e w \vi\sen{6 r) and vL = e w \vi\ cos(6 1 ^). Here e w , is a restitution 
coefficient characterizing the energy loss during collision with the walls, \vi\ 
is the magnitude of the pre-collision velocity and 6 R G [0, n] is a random 
angle @, M, EH • 



Equations of motion are integrated using a second order Velocity Verlet 
scheme [l2|, 13] and the results are given in dimensionless quantities using 
the normalization parameters m, g y d. Considering this, the unit for mass 
is m = 1, acceleration ^ = 1, diameter d = 1, time t = ^/d/g, angular 
frequency wq = \J g/d, velocity v = y/gd, force F = rug, elastic constant 
fco = mg/d, coefficient of viscous friction 70 = rriyjg/d and density po — 
m/d 2 . Table [21 lists this quantities and the values used in our simulations. 
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3. Results and discussion 

In all our simulations, the amplitude and frequency of the vertical accel- 
eration are a y = 6.83 d and u = 0.49 ujq respectively. As we mentioned in the 
introduction, the height difference Ay reaches a saturation value Ay sat after 
vertical vibration is applied. In Fig. [2] we shown the evolution the height 
difference Ay, where an accumulation of grains towards one side of the tube 
happens without any interstitial fluid, in a simple two dimensional setup and 
in the absence of an horizontal component in the oscillation. 




— Ay(t) = Ay sat [1 -e*-" 1 ' 1 ] 
Ay saf = 10.76, Tj = 1139.21 
J i i i i I i i i i l_ 
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Figure 3: Evolution of Ay for two different initial conditions, (a) Starting with the arms 
leveled, (b) Starting with an initial height difference of Ayo = 20,86<i. Segmented lines 
are saturating exponential fits discussed in the text. 



As we mentioned previously, the low accumulation of grains at reduced 
air pressure was initially observed in the experiments of Akiyama et al. [if. 
However, there is the reasonable doubt that the small accumulation observed 
by Akiyama et al. could be due to a small horizontal component present in 
the oscillation (since they do not discuss the issue of the verticality of the 
oscillation in their work). The presence of a small horizontal component can 
trigger the motion of grains to one side of a tube, as shown in reference (2)]. 
The simulations of Ohtsuki et al. |6|, in principle, solved this issue, because 
the oscillation is purely on the vertical direction, but their system is not sym- 
metric since their partitioned container has one branch wider than the other. 
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In our case, we observed the accumulation in a perfectly symmetric tube vi- 
brated in a pure vertical way This points out to the fact that the ingredient 
behind the accumulation is present in our simulations. It is also noticeable 
that the accumulation is observed in both the simulation of Ohtsuki et. al 
and ours, even with the different energy input mechanisms (Ohtsuki et al. 
shake the tube with a triangular waveform, while we shake used a sinusoid 
waveform) and different wall grain interaction (Ohtsuki et. al. simulated 
smooth walls, whereas we simulated rough walls). 

The saturation value of Ay, is quite robust with respect to different initial 
conditions, if we begin the simulation with the arms leveled (Ayo = 0), after 
a time t = 23400 1 Ay grows and reaches a saturation value of Ay sat rs 
lid (see Fig. 13^). If we start the simulations with the arms unbalanced 
(Ay = 20.86 d), we found that Ay decreases down to a saturating value of 
Ay sa t ~ 11 d (see Fig. [3b) that is the same saturating value for the growing 
case. 



14 1 — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




Figure 4: Saturation height difference Ay sa t as a function of the tube width dt, showing a 
peak for intermediate widths. Each point is the average of ten runs started with the arms 
leveled. The parameters characterizing grain-grain and grain-wall interactions were fixed 
to p = 1.27 po, k n = 20000 k , 7n = 20 7o , 7s = 20 7o , e w = 0.9. 

Our results for the evolution of Ay are well fitted by saturating expo- 
nential functions for both cases depicted in Fig. [3j The growing case (a) 
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is fitted by Ay(t) = Ay sat [l — e~*/ Tl ] with Ay sat representing the saturation 
value reached and T\ is the characteristic growth time. The decreasing case 
(b) is well fitted by Ay(t) = Ay sat [l + e - '/ 7 " 2 ]. The fact that T\ > r 2 can be 
understood since the transport of grains in the decreasing case is favored by 
gravity. The quantities djr\ ~ 10~ 4 and rf/r 2 ~ 10~ 3 , gives us an order of 
magnitude of the local effective velocity of a single grain moving in the tube, 
while Ay sat /T\ ~ 10~ 3 and Ay sat /r 2 ~ gives us an order of magnitude 
of the velocity of the bulk motion of the granular material through the tube. 
Both velocities are much lower than the velocity vo = y/gd = 1 of a single 
grain in free fall after covering a distance equal to its own diameter. 

The saturation height Ay sat is the main focus of our subsequent results, 
since it is a effective measurement of the amount of accumulated grains in 
one side of the tube. The fact that the accumulation takes place up to a 
maximum height difference Ay sat , that is the same whether the simulation is 
started with the arms of the tube leveled or with an initial height difference 
larger than Ay sat , implies that the growing process of Ay stops at the same 
point as the decreasing process, and further sustains the importance of the 
final height difference as a relevant parameter to monitor. 

14 
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Figure 5: Saturation height difference Ay sa t as a function of the grain density, for a fixed 
tube width of dt — lid, k n = 20000 fco, 7„ = 2O70, 7s = 2O70, e w = 0.9. Each point is 
the average of ten runs started with the arms leveled. 
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Figure S] shows the influence of the width of the tube d t on Ay sat . Each 
point is obtained from several simulations starting with the arms leveled. For 
each tube width, the number of grains iV was adjusted so as to have a fixed 
initial height of grains in each branch of L y /2. We can appreciate that the 
lower and larger tube widths (d t = 5.5 d and dt = 16.7 d) show practically no 
unbalance. This result suggest that there is an optimal tube width around 
dt = lid that maximize Ay sat . This result can be understood since for very 
narrow tubes, the flow of grains may be severely affected by the probability of 
an obstruction occurring inside the tube, thus no accumulation develops for 
very narrow tubes. The fact that there is no accumulation in the other limit 
(large widths), is understood since in that case, the walls through a complex 
mechanism linked to the vibration and the friction force, cannot maintain the 
weight of the columns, therefore is not produced the accumulation of grains. 




Figure 6: Ay sa t as a function of the grain- wall restitution coefficient e w . These observation 
indicates that under conditions of large energy dissipation, the accumulation is frustrated. 
This results agrees with previous results reported in [f| Q] . 

To probe the effect of the grain density on the accumulation, we plot in 
Fig. the value of Ay sat for a fixed tube width varying the grain density 
to change the weight of the columns. The simulations were done using a 
tube width of dt = 11 d, value around which the accumulative effect is more 
prominent, according to the results shown in figure HJ The value of Ay sat 
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decreases with increasing grain density. This result points out the fact that 
the accumulation process benefits from a low density granular material. At 
first glance, one is prompted to believe that the tube can hold the growing 
column only if its weight remains below a certain value. This behaviour is 
similar to an equilibrium condition, were the frictional force between wall 
and grains is not enough to hold up the granular column. 

To probe the effect of the wall on the growth of the columns, in Fig. 
[H] we show the variation of Ay sat with the restitution coefficient e w . Ay sat 
decreases rapidly when e w goes to zero. This effect can be understood as 
follows, the grains that collide with the walls lose kinetic energy abruptly due 
the low value of e w , these grains work like a "strongly dissipative granular 
trap" for the other grains, causing that lose energy quickly and stopping 
the accumulation process. For values of e w > 0.6, Ay sat reaches an almost 
constant value. In this case, the effect of walls on the grains is poor and the 
energy lost is due practically to grain-grain interaction, which is not changed 
during the simulations. 

4. Conclusion 

We studied the development of an accumulation of grains in one of the 
branches of a vertically vibrated two dimensional U-tube, in the absence of 
interstitial fluid and horizontal component. By studying the evolution of the 
height difference between the centers of mass of the arms of the tube Ay, we 
found that Ay reaches a saturation value Ay sat , that does not depend on the 
initial height difference. We probed the influence of the tube width, the grain 
density and the grain-wall coefficient of restitution on the development of the 
accumulation, by observing their effect on Ay sat . We found that there is an 
optimal tube width maximizing the development of the accumulation, and 
that it does not develop for too narrow or too wide tubes. The unbalance is 
promoted for low grain densities and frustrated at larger grain densities. If 
the coefficient of restitution is small, the accumulation is also frustrated, but 
if the restitution coefficient is increased, there is a point over which Ay sat 
becomes independent of it. 
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